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Abstract 

By making decomposition of the Wigner function simultaneously in both 

the spinor and the isospin spaces we derive a set of kinetic equations for the 

quark distribution functions and the spin densities. A detailed analysis of the 

consequences imposed by the chiral invariance on the form of the transport 

equations is presented. 
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I. INTRODUCTION 



The purpose of this paper is to generahze the already existing formulation of the mean- 
field transport theory for the Nambu - Jona-Lasinio (NJL) model The main aspect of 
this generalization is an extension of the simplified one-flavor approach to the more real- 
istic two-flavor case. Formally, this is achieved by applying the technique of decomposition 
of the Wigner function in both the spinor and the isospin spaces. In this way we can also 
extend and complement some of the earlier calculations done in the framework of QED |^ 
and QHD (quantum hadrondynamics) where the Wigner function was decomposed only 
in the spinor space. 

In the present approach we take into account all the coefficients in the spinor and isospin 
decomposition. Thus, although our formulation is restricted to the mean-field approximation 
it does not include any further simplifying assumptions concerning the structure of the 
Wigner function. This allows us to study the space-time evolution of the quark distribution 
functions and the dynamics of spin. We include also the possibility of having nonzero 
pseudoscalar condensates, which is crucial in studies of the chiral invariance of the theory. 

We concentrate in more details on the two special cases. The first one corresponds to 
the exact chiral limit, i.e., to the situation when the current quark masses vanish. On the 
other hand, the second case corresponds to the physical situation when the current quark 
masses do not vanish and are (slightly) different from each other. Studying the chiral limit we 
investigate how the concepts of chiral invariance can be explicitly included in the formulation 
of the transport theory. In this case we follow the treatment of Ref . 0] . 

The study of the transport theory for the NJL model becomes interesting in the context of 
the ultra-relativistic heavy-ion collisions. One expects that these highly energetic processes 
offer the possibility of creation of a short-lived quark-gluon plasma (QGP) 0. It is very 
likely that during such a deconfinement phase transition the chiral symmetry is additionally 
restored. The last fact is indicated by the lattice simulations of QCD ||^ showing that the 
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two phase transitions occur at the same temperature. Transport theory based on the NJL 
model gives us the possibihty to study phenomena connected with the chiral phase transition, 
which take place in the systems out of thermodynamic and chemical equilibrium. This is 
an attractive feature of the model, since the methods for applying QCD directly to such 
situations have not been elaborated yet. The QCD transport theory has been formulated 
in papers by Heinz |§] and by Elze et al. P| but this approach is not frequent in practical 
applications. 

In the mean-field approximation the NJL model includes only the quark degrees of free- 
dom. Therefore, this approach is not fully appropriate for the description of hadronic matter 
at low temperatures or densities and one has to go beyond the mean-field approximation in 
order to obtain the agreement of the low-temperature NJL results with chiral perturbation 
theory |10|. On the other hand, close to the deconfinement and chiral phase transitions 
the quark degrees of freedom become relevant. In this situation the mean- field approach is 
appropriate since it describes many important features of the chiral phase transition like, 
e.g., the decrease of the in-medium quark condensate. Consequently, the present formulation 
of the transport theory is most suitable for description of the phenomena happening in the 
neighborhood of the chiral phase transition. Work concerning inclusion of the meson degrees 
of freedom into the transport theory for the NJL model (this is equivalent to the extension 
of the mean-field approach) is being currently carried out by the Heidelberg group 



The paper is organized as follows. In the next Section we define the model. In Sect. |T1 
we introduce the Wigner function, define its spinor and isospin decomposition, and discuss 
the chiral transformation rules. Sect. |^ presents the quantum kinetic equations satisfied 
by the coefficients of the spinor and isospin decomposition. In Sect. we do the classical 
approximation and derive the so-called constraint equations. We do it separately for massless 
and massive quarks. The classical kinetic equations for the quark distribution functions and 
for spin densities are derived in Sect. [VJ. We summarize in Sect. |V11| . 
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II. DEFINITION OF THE MODEL 



In this paper we take into consideration Lagrangian 

G 
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{^>mY + {^>i^^T^Y , (1) 



where \1/ = {ipu-, V'd) is the doublet of the Dirac fields, G is the coupling constant, 
T = (Ti,r2,T3) are the Pauli isospin matrices, and rh is the matrix containing current 
quark masses ttIu and md. This matrix can be written in the form rh = rriQ + m-r, where 
rriQ = ^{rriu + rrid), mi = 1112 = and = |(m„ — m^). For simplicity, the color degrees of 
freedom of quarks have been neglected. 

Using Lagrangian (|l|) and assuming mean-field approximation we find the following equa- 
tion for the field 

[i^ - m - a{x) - i'j57r{x)-T]'^{x) = 0. (2) 
Here the mean fields cr(x) and 7r(x) are defined through expressions 

a{x) = -G ( ^{x)^{x)) = -G Tr (p ^{x)'^{x)) (3) 

and 

7r(x) = -G ( ^{x)i^5T^{x)) = -G Tr (p ^(a;)i75r^(x)) , (4) 

where p is the density operator and Tr denotes the trace over the Fock space. 

The important feature of Lagrangian (|l]) are different symmetries which resemble us the 
case of QCD. First of all Lagrangian (|l]) is invariant under Uv(l) transformations, which 
leads to the conservation of the baryon current 

df.Vf'ix) = 0, V^'ix) = {^{x)Y'^{x)). (5) 

In the isospin symmetric case, m„ = rrid, Lagrangian (|ip is additionally invariant under 
SUv/(2) transformations. This fact gives the conservation of the isospin current 



d^yix) = 0, Y^{x) = (^(x)7''t^(x)). (6) 

In the chiral limit, m„ = = 0, expression (|l|) does not change under SU^(2) transforma- 
tions. This invariance leads to the conservation of the axial current 

d^A'^ix) = 0, A'^(x) = (^'(x)7^75T^(x)). (7) 

We note that we have defined the conserved currents as the expectation values in the state 
characterized by the density matrix p. In this case the conservation laws follow directly from 
Eqs. (i - ®. 

Let us now introduce the Green function 

where i,j and a, (3 are isospin and spinor indices, respectively (i, j = 1,2 and a, /3 = 1, 4). 
One can easily check that the Green function (|]) fulfills the same equation as the field \& 
does, namely 

[i ^ — m — a{x) — i757r(x) -r] S^{x, y) = 0. (9) 

Moreover, the mean fields a{x) and 7r(x) can be determined directly from S^{x,y) via 
relations 

a{x) = —G ti S'^{x,x), 7r(x) = — G tr i75TS'^(x, x), (10) 

with the trace tr taken over spinor and flavor indices. One can observe that formulas and 
(p!0|) form a closed system of equations. It will be the subject of our studies in the following 
Chapters. 

III. WIGNER FUNCTION 

A. Decomposition in spinor and isospin space 

Our aim is to derive a set of kinetic equations for the classical distribution functions. 
This can be achieved in the usual way by introducing the Wigner function 
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w. 



{X,p) = ^J d\ei^-S<^^^{x + lu,X-lu) (11) 



ija/3 

and by searching for equations fulfilled by this function in the limit h —>■ 0. In Eq. ( pT]) the 
quantity X is the center-of-mass coordinate, X = ^{x + y), and u is the relative coordinate, 
u = X — y. 

Starting from Eq. and using the well known results for the Wigner transform of the 
derivative of a two-point function, df{x,y)/dx^, and the Wigner transform of the product 
of a one-point function with the two-point one, f{x)g{x,y), we find the following equation 



- M(X) - m-T + '4^^^M{X) - 275 7z{X)-t - ^75 d^7ziX)-Td^ 



W{X,p) = 0, 
(12) 



where we used the standard notation = d/dX^, = d/dp^, = + and 
M{X) = (t{X) + niQ. We note that in Eq. (0) higher order gradients have been neglected, 
so it is valid only for weakly inhomogeneous systems. In addition, there is no collision term 
on the right-hand side of Eq. ([I2|) , which is a consequence of our mean- field approximation. 



The effects of collisions have been recently studied in Ref. [|I2|. 

Since the Wigner function satisfies the condition W{X,p) = 7°iy^(X, p)7° = W{X,p), 
it can be represented as the following combination of the Dirac tensors F 

W = :F + i^r^V + i^V, + n^A, + \a'''S,,. (13) 

In the decomposition (p!3|) each coefficient C [i.e., the functions JF(X, p),'P(X, p),V^(X, p), 
A^{X,p) and iS^jy(X, p)] is a hermitian two by two matrix. Thus, it can be further decom- 
posed in the isospin space according to the rule 

C = C + C-T. (14) 

In this way we generalize the usual approach, which does not take into account such 
structure. Similarly to Eq. (|I^) we introduce the quantity M{X) defined as M{X) = 
M{X) + m-T = cr{x) + rriQ + m-r. 
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Many of the functions defined by the spinor and isospin decomposition, Eqs. (|T^) and 
(0), have a direct physical interpretation. For example, the quantities V'^{X,p),V^{X,p) 
and ^^{X,p) are the phase-space densities of the baryon, isospin and axial currents. More- 
over, according to Eq. (p!OD, the mean fields (^(X) and 7r(X) are simply related to the 
functions jF(X,p) and 7^{X,p) 

a{X) = -8Gj^^:F{X,p) (15) 



and 



7r{X)=8Gj-0^^nX,p). (16) 
For the physical interpretation of the other components we refer the reader to Refs. p|,[T3 . 



B. Chiral transformations 

The most important feature of Lagrangian (|I]) is its exact chiral invariance in the case 
m = 0. Studying the transport theory it is interesting to analyze the consequences of this 
symmetry for the transport phenomena. This has been already discussed in detail for the 
one-flavor version of the model P]. In this paper we shall generalize these results to the 
two-flavor case. 

The SU^(2) chiral transformation of the fleld \E' is deflned as follows 

^ ^ ^' = exp(-i75^^)^ = ^cos|^ -i75n-Tsin0 ^. (17) 

Here n is a unit vector in direction of x ^iid x is the length of X- The property (|l^) induces 
the transformation rules: 

(T — >■ o"' = (T cos X — Ti" ■ n sin X (18) 

for the scalar mean fleld and 

TT — i> tt' = TT cos^ ^ — [2(7r-n)n — tt] sin^ ^ -|- an sin x (19) 



for the pseudoscalar mean field, respectively. As is expected, Eqs. (|T8|) and (|T9D indicate 
that the sum of squares o"^ + tt^ is an invariant of the chiral transformations. Another 
consequence of Eq. (0) is the transformation law for the Wigner function 

W ^W' = exp(-^75^) W exp(-275^). (20) 

Eq. ( pOD leads to a set of the transformation rules for the coefficients in the spinor and isospin 
decomposition. Their full form is listed in the Appendix. Below we give the prescriptions 
for the infinitesimal chiral transformations (denoting the infinitesimal value of x by 5%). 

Infinitesimal chiral transformations of the scalar and pseudoscalar coefficients: 

jr^jr' = jr + p.n^;^, (21a) 

^ J"' = T + Vn5x, (21b) 

V -^V' = V - T-n5x, (21c) 

r ^V' = V -J^n6x. (21d) 

Infinitesimal chiral transformations of the vector and axial-vector coefficients: 

(22a) 
(22b) 
(22c) 
(22d) 

The quantities S^i,,S^u,Sfj,i, and S^u [the last two are the dual spin tensors defined below 
in Eq. (plj) ] transform in the same way as the functions JF, V and 7-*. This property 
resembles the case of the one-fiavor calculation (see Eq. (22) of Ref. [Q). In addition, 
similarly to the one-fiavor case one finds that both and are chirally invariant. 
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IV. KINETIC EQUATIONS 



Substituting formula ( P^ into Eq. (^) and comparing the coefficients appearing at the 
Dirac tensors we find a set of coupled equations 
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(23a) 



(23b) 



(23c) 
(23d) 



(23e) 



In Eqs. ( f^3d| ) and (|23e| ) we have introduced the dual spin tensor S defined through relation 

S =\e^^^^S^p. (24) 

Equations (|23| ) represent the spinor decomposition of Eq. ([T^). They are a generalization 
of the set of equations (24) - (28) from Ref. ||^ where a one-fiavor formulation has been 
analyzed. In our two-fiavor approach each of Eqs. ( p3D has a matrix form. In order to 
obtain a system of equations that is easier to control, we perform further the isospin de- 
composition. This procedure allows us to deal only with real quantities, which represent 
physical observables. At the first stage of the isospin decomposition we insert expressions 
of the form ([I^ ) into Eqs. (|23D . Subsequently, we calculate the sum and the difference of 



the initial matrix equation and its adjoint for all the formulas appearing in (|23|) . In the 
resulting five pairs of equations one compares the coefficients at the Pauli isospin matrices. 
In doing so, it is convenient to use two relations 

z-{t,C} =2z-C + 2Cz-T (25) 

and 



z-[T,C]^2i{zxC)r, (26) 

where z is an arbitrary three- vector and [,]({, }) denotes the commutator (anticommutator). 
Finally, we get five groups of equations. 

Scalar equations: 

p^'Vi.-MJ^-m-J^+n-V^O, (27a) 



-'^V^-M:F-mJ^ + 7rP = --(a,7rxa;p), (27b) 



h 

p^V^ - M - mj- +7rV = — 

^d'^v, = -\d.Md;r + |a.7r-a;p, (27c) 

^a'^V^ - mx:F + 7rxP = -^a.M^^JF + ^d,Tzd;V. (27d) 



Pseudoscalar equations: 

n 



d^A,,- MV -m-V -TT-J^^O, (28a) 

\d^A^ -MV-raV-T^r^ ^d.n x d;j^, (28b) 

p^A = I^.Afa;^' + |9,7r-a;:F, (28c) 

inxV + nxJ^ + jfA^ = ^9,M9;p + ^9,7r9;jF. (28d) 

Vector equations: 

p^T - h^S,^ - MV^ - m- = ~d,7T-d;A^, (29a) 

-nxA^+Pi,J^- \&'S,^ - M V;. - mV^ = -^^.Tr^;^^, (29b) 

^d^J=^p-S,^ + TT-^^ = _^a.Ma;V^, (29c) 

- m X V;. + \d^j= + + TT^^ = -^a,Ma;v^ - ^a,7r x a;^^. (29d) 

Axial- vector equations: 

- \d^V - pj""^ - MA" -nv A"^ ~d,n-d;V", (30a) 
- TT X - ^d^V - PuS""" - MA^" - m^^ = -^9,7r-5;V^, (30b) 

p'^V - hj'"' + 7r-V^^ -^-d,Md;A", (30c) 



m 



2" 2 

X + p^'v - + ttV^ = ~d,M d;A" - \d,Tz X a; v^ (30d) 

10 



Tensor equations: 



(31a) 



(31b) 



(31c) 



(31d) 

The first two equations in each set of formulas ( ^Tf ) - (0) can be regarded as the "her- 
mitian" parts of Eqs. (p3D. They correspond to Eqs. (31) - (35) of the one-flavor approach 
0]. The other equations in formulas ( pTf ) - ( |3lD are the "anti-hermitian" parts of Eqs. ( p3D 
and correspond to Eqs. (37) - (41) of 0. In the case m = 0, using expressions (^) and 
(p^) defining infinitesimal chiral transformations, one can check that Eqs. ( pTf ) - ( pT]) are 
chirally invariant. 

One can notice that Eq. ( ^7c| ), integrated over four-momentum p, leads to the baryon 
current conservation (|^). In the analogous way, using Eqs. (|27d|) and ([T6|) one finds the 
formula 

hd,\^{X) = 16m X I 0-^J^{X,p), (32) 

which yields the conservation of the isospin current for the symmetric case = [compare 
Eq. (P)]. Finally, Eqs. (|28b| ), (|^) and (|16D lead to the expression 

hd.A^iX) = -^7r(X) - 16m| (33) 

which is reduced to the axial conservation law, Eq. (|^, in the chiral limit m„ = rrid = 0. 
Thus we see that after making the gradient expansion, the conservation laws are still included 
in the transport equations. 
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V. CONSTRAINT EQUATIONS 



A. Classical Approximation 

In order to obtain classical transport equations one makes an expansion of the functions 
^{X,p), V{X,p), V^(X,p), Afj,{X^p), S^u{X,p), <7{X) and 7r(X) in powers of h. In this 
way each coefficient in the decomposition (|T3|) as well as the functions a and tt can be 
represented as a series 

C = C^o) + hC^i)+n'C(^2)+ ... . (34) 



Inserting expressions of the form ( p4|) into the kinetic equations (^) - (0) and comparing 
the terms appearing in the leading (zeroth) order of h we find a set of constraint equations, 
which connect different leading order terms of the coefficients C. 

Scalar constraint equations: 

pf^yio) _ M(o)^(o) - m-:F(o) + 77(0) •P(o) = 0, (35a) 

Vjf) - M(o) :^(o) - m ^(o) + tt (o) P(o) = 0, (35b) 

- m X :F(o) + 77(0) X V(o} = 0. (35c) 

Pseudoscalar constraint equations: 

M(o)V(o) + m-P(o) + 77(0) -^(0) = 0, (36a) 

M(o) P(o) + m P(o) + 77(0) ^(0) = 0, (36b) 

= 0, (36c) 

m X P(o) + 77(0) X :F(o) + p'^^Jf ) = 0. (36d) 

Vector constraint equations: 

p^^(o) - M(o)V(°) - m- V(°) = 0, (37a) 

- 77(0) X + p, :^(o) - M(o) V(°) - m = 0, (37b) 
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p^5(°)+7r(o)-4°^ = 0, 
-mx Vf+p^5W + 77(0)4°) = 0. 



(37c) 
(37d) 



Axial-vector constraint equations: 



pJl^^ + M(o)4o) + ^--^U = 0' (38a) 

^(0) X Vfo) + + M(o)^fo) + m 4) = 0' (38b) 

P^7^(o) + 7r(o)-Vfo) = 0, (38c) 

- m X Afo) + p'^P(o) + 77(0) Vfo) = 0. (38d) 

Tensor constraint equations: 

- e^^'^-^p^Af - 77(0) -^J^) + M(o)5(';^ + m-5;5 = 0, (39a) 

- e^^-^p^Af - 77(0)4';'; + ^(O-^W + ™ '^S) = 0' (39b) 

p^Vl,) - p^Vf^) = 0, (39c) 

- 77(0) X 5f;) + m X S^^^ + p^'V^o) - P^Vfo) = 0. (39d) 

The constraint equations written above have to be supplemented by the formulas which 
determine the values of the mean fields in the leading order of h 

d'^p 



a(o) (X) = -8G / ^^(0) {X, p) (40) 



and 



77(o)(X)=8g/-^P(o)(X,p). (41) 



(27r) 



B. Chiral limit 

In this subsection we shall analyze the form of the constraint equations ( pSj ) - (^) in 
the limit itlq = m = [in this case M(o)(X) = o"(o)(X)]. Having in mind the previous 
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studies relying on the spinor decomposition we expect that quantities ^(o)) -^(o)) 

A'^Q'j and ^^Q-j can be used as the fundamental variables in construction of the transport 
theory: JF(q) describes the quark space-time distribution, ^(q) specifies the isospin of quarks, 
Al^Q-j describes the quark spin density, and At^^-^ determines the spin density of quarks with 
different isospin components. We note that the "classical" isospin of quarks is described here 



by a three-vector. This fact is connected with the form of the decomposition (|Tj) of the 
two by two, quantum mechanical density matrix. We also note that Eqs. (|36c|) and ( p6d|) 
indicate that only three out of four Lorentz components of A'^q^ and A!^q^ are independent. 
Therefore, similarly to the "classical" isospin, the "classical" spin of quarks is also described 
by a three-vector. 



Let us now express the quantities P(o), 'P(o), V^g^, Vj^g)' '^(o)' '^(o)' •-'(o) 



S'f'r'^-, and S^Q-j as 



(0)) -'"(0); -^(0) ^iid ^['q). First of all, it is easy to notice that Eqs. (|36aj) and 



functions of J-'{(^^, J-^ 
( p6b|) define the pseudoscalar densities in terms of T{q) and J- (q) 

7r(o)-^(o) 



(0) 



(42) 



'(0) 



(O)-^(O) 



(43) 



In the similar way Eqs. ( P7a| ) and ( |37b|) define V^^-, and V('q) in terms of T(^) , ^(o) and M^^^ 



•^(0) - V 



^(0) 



(44) 



»^(0) - V 



(0) 



77(0) X A 



•(0) 



^(0) 



(45) 



After a few algebraic manipulations Eqs. ( p9a| ) and ( |39b| ) allow us to express the spin tensors 



S'^Q^ and the dual spin tensors Sf^^^, S(q) as functions of A'^q) and Al^ 



cr 



(0) ^iLvafi 



PaA 



(0) 
/3 ) 



(46) 



P ^(0) P -^(0) 



A I 



(0) 



(47) 
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(0) 



P ^(0) P '^(O) ^ ^ P<^ 



(0) 



(0) 



(0) 



(48) 



and 



TT 



(0) 



p 



TT 



cr, 



(0) 



(0) 



(0) 



■/3 



o-(o) 



(0) P -^(O) 



(49) 



Here we have introduced the chirally invariant mass 



^2(X)=a(2o)(X) + 7rfo)(X) 



(0)> 



(50) 



Let us now take into consideration other constraint equations appearing in the leading 
order of h. Substituting expressions ( ^2]) - (^Sf ) into Eqs. (|35a|) and (|35b| ) one finds two 
mass-shell conditions 



[p'-M'iX)]J^^o)iX,p) = 



(51) 



and 



[p'-A^2(x)]:F(o)(x,p) = o. 



(52) 



The third scalar constraint equation is automatically fulfilled, since 7^(0) and 7r(o) are parallel 



see Eqs. ( p5c|) and (H)]. 



Substituting expression (|^) in Eq. ( p7c| ) and using condition (|36d| ) one finds 

[p'-M'{X)]n^oyA1,^=0. 
In the analogous way, substituting (|^) in (|37d| ) and using condition ( P6c|) we get 



(53) 



[p'-M'iX)]A1,^ = 



(54) 



This formula follows also from Eqs. ( p8a| ), (^) and ( |36b| ). On the other hand, Eqs. ( |38bD , 
(H3), (BB3) and (E3|) lead to the mass-shell constraint 
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[p'-M\X)]A^,,=0. (55) 



•-(0) 

To complete our discussion of the constraint equations in the chiral hmit we should 
discuss an important property of formulas (^) and (^). Using Eq. ( ^31) in (^Tj) we find 
that Eqs. (^) and (^) are not independent — they do not determine separately the values 
of (7(0) (X) and 7r(o)(X). This fact is a consequence of the chiral symmetry. It indicates that 
only the invariant mass Ai{X) has a physical significance. 

C. Massive quarks 

In this subsection we are going to consider the case m 7^ 0. Calculating 7^(0) from Eq. 
( p6b|) and substituting into Eq. (^) we find that 7r(o) and m must be parallel in this case. 



Using again Eq. ( p6b|) we find that 7-*(o) must be also parallel to m. Further inspection of 



Eqs. (|35|) - (|39|) indicates that all quantities C(o) are parallel to m [i.e., only their third 
component is different from zero, C(o) = (0,0,C(o)3)] . In this situation it is convenient to 
introduce the combinations 

C(o)u = C'(o) + C'(o)3, C'(o)d = C(o) - (^(0)3, (56) 
which describe the up and down quarks, respectively. Using this notation, we can write 

no) u = -7r(o) 3 , . , , _ , no) d = 7r(o) 3 — ::- (57) 



M(o) + m3' ^"^^-^ M(o)-m3 



and 



''W'^-M(o) + m3' ''(°)'^-M(o)-m3- ^'^^ 
In the massive case the gap equations (|40|) and (^Tj) take the form 

J (27r)4 1^ a(o)(X) a(o)(X) ) ^ > 



and 
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7r(o)3(^) 



1 + 46- 



0. 



(60) 



(27r)4 ycr(o)(X) + rriu cr(o)(X) + m^^ 
These two equations lead us to the condition 7r(o)3(X) = 0. Thus, in the case when the 
chiral symmetry is exphcitly broken, the leading contribution to the pseudoscalar condensate 
should vanish (similarly as in the one- flavor case, see Eq. (86) of 0]). 

Substituting expressions (|58D into the scalar constraint equations ( P5a| ) and (|35b| ) , we 
find the mass-shell conditions 



- (M^(X))2] ^(o)/(X,p) =0 (/ = u,d), 
where we have introduced the notation 

M"(X) = M(o)(X) + ma, M'^(X) = M(o)(X) - ms- 

We turn now to the discussion of the spin dynamics. First of all, Eqs. ( |36c| ) and 
give us the condition 



(61) 



(62) 



(0)/ 



0. 



(63) 



The tensor constraint equations ( |39aD and ( |39b| ) can be used to find the formulas for the 
spin tensors, namely 



(0)/ 



(64) 



The dual spin tensors are obtained by contracting this expression with the Levi-Civita tensor 
as in (0). Using now Eqs. (je|), (|B|), (pgiD and (§S§) we find 

'p'-{Mf{X)r]A^,^f{X,p) = 0. (65) 



VI. CLASSICAL TRANSPORT EQUATIONS 

In this Section we derive classical transport equations. This requires the study of Eqs. 
(p7|) - (0) up to the first order in h. Since the calculations done in the chiral limit have 
different aspects from those done for massive current quarks, we discuss these two cases 
separately. 
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A. Chiral limit 



Substituting Eqs. (||) and (|§ into Eq. we find 



p'^df.F + Md^Md^F = 0, 



(66) 



where F{X,p) = jF(o) (X, p) /cr(o) (X) is the chirally invariant quark distribution function 
(compare Eq. (60) of 0). In the similar way, substituting (^2]), ( ^3D and the formula for 
'P{X,p) obtained from ( |28b| ) into Eq. ( ^7d| ) one gets 



(0) 



cr. 



(0) 



(67) 



where F{X,p) = ^(o)(X, p)/(T(o)(X). In contrast to F{X,p) the function F{X,p) is 
not chirally invariant [under infinitesimal chiral transformations F — > F' = F + 7r(o) x 
(F X n) (^(5x/o'(o))]. However, one can check that the form of Eq. ( |67D is chirally invariant. 

In order to find the kinetic equations satisfied by the functions A^^^-^ (X, p) and ^^^^ (X, p) 
we use Eqs. (|30c|) and ( |30d| ), respectively. Calculating V{X, p) from Eq. ( P8a| ) and V^'q) (X, p) 
from Eq. ( p9b| ), and substituting these two expressions into Eq. ( pOc| ) we find (to the first 
order in h) 



(68) 



In the analogous way, calculating 7-*(X, p) from Eq. ( p8b|) and V^(X, p) from Eq. (p9aD 
and inserting these expressions into Eq. (|30d|) we find (again to the first order in h) 



(69) 



Using now expressions ( ^61) - (^91) , defining the leading order parts of the spin tensors, we 
find the desired equations 



p-^d^A^,^ + Md^Md^A^,^ + — 



P >l(Q) P ^(Q) 



-e^-^p^—dAm\-J^'=Q 



M ) '^^ 



(70) 



and 



P -^{0) P -^(0) 





■/3 - 



to) -P'''^{0)-^(C 



(7(0)9^77(0) X F. 



+ 



(71) 



Several comments are in order now: 



(a) One can check that Eqs. (|70D and (|7TD are consistent with the conditions p'^A^^^ = 



and p^A^^' = ^(0) X 7r(o). One can also check by a straightforward but tedious calculation 
that Eqs. (|70D and (^) are chirally invariant. 

(b) The kinetic equations (|66D , (|67|), (|70D and (|7TD together with the mass-shell con- 
straints (|5T|) , (|5^) , (|5^) and (^) describe the evolution of quarks in the given scalar and 
pseudoscalar fields. The calculations within the NJL model require that these fields are 
obtained self-consistently from the quark distribution functions. However, the discussed set 
of equations can be also used in the cases where we are interested in the dynamics of quarks 
in the externally given scalar and pseudoscalar fields. For that reason, the results presented 
in this Section are an extension of a part of the earlier results obtained in |l4| , p!5| . An in- 
teresting novel feature of the present approach is the coupling between the spin and isospin 
degrees of freedom. 

(c) Doing calculations in the NJL model, one evaluates the mean fields from the self- 
consistent equations (|l^) and (|T^) . However, due to the chiral symmetry of the model [see 
discussion following Eqs. (|55|)] only a chirally invariant combination Ai^ = 7r(o) + o"(^o) "-^^^ 
calculated from the knowledge of the quark distribution function. Thus, the discussed system 
of equations is closed (and can be solved) only in the case when the ratio 7r(o)(X)/cr(o)(X) 
is fixed and independent of the space-time position coordinate X. Nevertheless, this ratio 
can be arbitrary, which is required by the chiral invariance of the system. 

(d) Since the ratio 7r(o)(X)/cr(o)(X) is fixed and arbitrary, in the practical calculations 
we can always set it to zero (this is like choosing the convenient gauge in the QED-type 
calculations). In this case the form of the kinetic equations (^), (|70| ) and (|7l| ) simphfies 
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substantially and reduces to that known from the one-flavor approach (compare Eqs. (62) 
and (68) of Ref. [|). 

(e) Eq. (|7I]) does not guarantee the conservation of the axial current, see Eqs. (|^ 
and (^). This fact can be seen in the special case, 7r(o)(X)/(j(o)(X) = 0, by using similar 
arguments to those presented in the one flavor study ||^. Consequently, solutions of Eq. 
([n]) are constrained by the condition of the axial current conservation. Only those solutions 
which satisfy Eq. (|^) are valid. 



B. Massive quarks 

Using Eq. ( p8b|) and presenting similar arguments to those used in Section 0, we find 
that the cross product tt x T' vanishes up to the first order in ft. In this situation Eqs. ( |27c|) 
and ( |27d| ) lead to the kinetic equations satisfied by the distribution functions of the up and 
down quarks. They have the following form 

Tfd^Ff{X,p) + Mf{X)d^Mf{X)d^Ff{X,p) = (/ = u,d), (72) 

where F^(X,p) = J^(o)/(X,p)/M^(X). 



The kinetic equations for the spin densities follow from Eqs. ( |30c|) and (|30d| ). Using the 
formulas for the pseudoscalar densities V and "Ps, calculated from Eqs. (|28a]) and (|28b|) up 
to the first order in Ti, we can rewrite ( pOcj ) and (|30d| ) in the form 



Ip^^d'^Af - ms (p'^P(i) 3 + vr(i) aVfo)) + \M^o)duSl^) = -\M^,)d,M^,)d;A^^,) (73) 

and 



Eqs. (p8a|) and (|28b|) give us additionally two conditions 



M«P(i) „ = i^.^o) u - ^(1) 3-^(0) u (75) 
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and 



(76) 



Substituting Eqs. (|58D, ([75| ) and ( |76| ) into the sum of Eqs. (|7|) and ([7^, and using formulas 
for the dual spin tensors we can find 



(77) 



In the similar way, substituting Eqs. (|58D , (f75|) and ([76|) into the difference of Eqs. (^) and 



(ff^) one finds 



(0), 



0. 



(78) 



Equations (^) and (|78D are the analogs of the spin kinetic equation derived for the first 
time in [Q]. 

The mass-shell conditions allow us to express F"''^(X, p) as the sum of the quark and 
antiquark distribution functions /J"^(X, p) and /~^(X, p) 



2E;'\X) 



2E;'\X) 



/„.(X,-p)-l (79) 



where 



e;^\x) = J(M«''^(x))2 + p2. 



Substituting this formula into Eq. (^), and integrating over gives 



p^9,/„±,(X, p) + M"''^(X)9,M"''^(X)9;/;,(X, p) = 0. 



(80) 



(81) 



Using expression (79) we can also rewrite the gap equation in the more familiar form 

d^p M"" 



M(o) - mo = -2G 
-2G 



(277)3 E«(X) 

d^p M'^ 
(2vr)3 E^(X) 



/+(X,p) + /„-(X,p)-l 
/+(X,p) + /,-(X,p)-l 



(82) 
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The last results indicate that the form of the kinetic equations for the two- flavor approach 
(with massive current quarks) is analogous to that of the one-flavour model. The up and 
down quarks, as well as their spins, evolve in the mean fields and M'^, respectively. This 
allows to introduce the spin up and spin down densities as the appropriate combinations of 
the functions {X,p) and A^Qy{X,p) (compare Eqs. 69 and 70 of 0). The only interplay 
between the two flavors occurs via the gap equation, which determines the common part of 
the mean fields M" and M"^. 

VII. SUMMARY 

In this paper we have derived and analyzed the mean-field transport equations for the 
two- flavor NJL model. In this way we have extended the previous approach 0, restricted 
to only one flavor. Applying technique of the decomposition of the Wigner function in both 
the spinor and isospin space, we could investigate the most general form of the mean-field 
classical transport equations. 

We have discussed in detail the case of a chiral limit, studying the invariance properties 
of the transport equations. Our analysis shows the limitations in the applications of the 
chirally invariant kinetic equations. This sort of restrictions is already known from the one- 
flavor considerations: the ratio 7r(o)(X)/(j(o)(X) must be constant and the spin dynamics is 
constrained by the axial current conservation. 

An additional outcome of our approach are the transport equations for quark matter 
moving in the externally given scalar and pseudoscalar fields. They follow from our analysis 
in the case when the self-consistency required in the evaluation of the meson mean fields is 
relaxed. These equations describe non-trivial couplings between the spin and isospin degrees 
of freedom. 

If the current quarks are massive, the difficulties connected with the requirement of 
the chiral invariance are not present. The space-time evolution of the up and down quark 
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distribution functions is determined by the mean fields M"(X) = M(o)(X) + |(m^j — m^) 
and M'^(X) = M(o)(X) — |(m„ — m^), respectively. The common part of the mean fields, 
M(o)(X), is obtained from the self-consistent gap equation. 
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APPENDIX: 

In this Appendix we list the chiral transformations rules obeyed by the coefficients ap- 
pearing in the spinor and isospin decomposition of the Wigner function (pT^). 



Scalar and pseudoscalar functions: 

^ JF' = jFcosx + "P-nsinx, (Ala) 

J= ^J=' = :Fcos^ I - [2(:F.n)n - J=\ sin^ | + Pnsinx, (Alb) 

p^p' = Pcosx-:?^-nsinx, (Ale) 

V = Pcos^ I - [2 (P -11)11 - P] sin^ I - J^^nsinx- (Aid) 

Vector and axial-vector functions: 

^ V; = V^, (A2a) 

^ V; = V^cos^ I + [2(V^-n)n - V^] - n x A^sinx, (A2b) 

^ = A^, (A2c) 

A^^ = cos^ I + [2(^^ ■ n)n - A^] - n X sin x, (A2d) 
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Tensor functions: 



S'nv = ^i^v COS X + ^^i. • n sin x, 

2 X 



s: 



= ^1^1^ COS x-S^y-n sin x, 



S^y COS- - - 



2(5^^-n)n - S^^ 



(A3a) 
(A3b) 
(A3c) 
(A3d) 
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